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Abstract 

We consider a single copy of a quantum particle moving in a potential and show that it is 
possible to monitor its complete wave function by only continuously measuring its position. While 
we assume that the potential is known, no information is available about its state initially. In order 
to monitor the wave function, an estimate of the wave function is propagated due to the influence 
of the potential and continuously updated according to the results of the position measurement. 
We demonstrate by numerical simulations that the estimation reaches arbitrary values of accuracy 
below 100% within a flnite time period for the potentials we study. In this way our method grants, 
a certain time after the beginning of the measurement, an accurate real-time record of the state 
evolution including the influence of the continuous measurement. Moreover, it is robust against 
sudden perturbations of the system as for example random momentum kicks from environmental 
particles, provided they occur not too frequently. 
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Monitoring - continuous observation - of a dynamical system in the presence of random- 
ness is employed not only in physics and chemistry, e.g., to survey the motion of comets, the 
growth of thin layers or the dynamics of chemical reactions. It is a sub-discipline of robotics 
and also plays a vital role in other fields such as earth sciences and aeronautics with numer- 
ous applications from climate observation, control of robots and vehicles to remote sensing. 
Monitoring tasks can be modeled by stochastic processes supported by continuous updates 
of estimates according to the observed random data - also called stochastic filtering A 
special challenge, however, is posed in the realm of quantum physics; the preparation - not 
to mention monitoring or control - of individual atoms, electrons and photons remained 
experimentally unattainable for half a century. The theory of quantum monitoring only 
emerged 20 years ago . Nowadays, however, nano-technology and quantum infor- 

mation processing strongly inspire a mathematical theory of monitoring and control of single 
quantum degrees of freedom like, e.g., the position of an atom or a nano-object. A more 
challenging aim is to monitor and control the entire state of individual quantum systems. 

Monitoring a quantum system encounters principle difficulties that lie in the characteristic 
traits of quantum nature itself: incompatible observables such as position and momentum 
and irreversible state change introduced by measurements. Methods have been developed 

n 

to employ monitoring in order to determine the pre-measurement state [6|, for parameter 



estimation 



purposes jll| or to reach a targeted state Il2[|. Moreover, the possibility of state monitoring 
has been studied for special systems 13|, . We are here going to show, that monitoring 
the position of a single quantum particle promises - via our theory - the monitoring of 
the full wave function, i.e., the complete state of a particle, and demonstrate its power by 
numerical simulations. Like always in the quantum realm, monitoring will unavoidably alter 
the original (unmonitored) evolution of the wave function. Strong monitoring assures very 
robust fidelity of the estimated wave function but it has little to do with the unmonitored 
wave function. Fortunately, in many cases, a suitable strength of monitoring assures both 
robust fidelity and slight change of self-dynamics. Needless to say, that such compromise is 
not due to any weakness in our theory. It is definitely an ultimate necessity enforced by the 
Heisenberg uncertainty relations. 

In the following the concept of continuous observation is interpreted as the asymptotic 
limit of dense sequences of unsharp position measurements on a single quantum particle. 



7|, to track Rabi oscillations 8|, [9| , Il0| and -combined with feedback- for cooling 
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We describe the inference of its wave function from the sequence of the measured position 
data and eventually compare the true wave function with its estimate using simulations of 
quantum particles moving in various potentials. Among them is the Henon-Heiles potential, 
which in classical physics implies chaotic behaviour and thus exposes tracking of dynamics 
to extreme conditions. 



I. MONITORING THE POSITION 

Time-continuous position measurement can be understood as an idealization of a se- 
quence of discrete unsharp position measurements carried out consecutively on a single copy 
of a quantum particle The notion of unsharp measurement is instrumental here. Such 
an unsharp measurement of the position q can be realized as indirect von Neumann mea- 
surement; instead of measuring the particle's position directly, an ancilla system is scattered 

17, Q, Q. The observed results 
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off the particle and then the ancilla is measured |15L 
yield limited information on the position q of the scatterer. In a simple description, a single 
unsharp measurement of resolution a collapses the wave function onto a neighborhood with 
characteristic extent a of a random value q: 

m — -^VGM^m , (1) 



where G„{q) = (1/^2^) exp(— g^/2(T^) is a central Gaussian function. The random quan- 
tity q is the measured position which determines the collapse, i.e., the weighted projection, 
of the wave function. The probability to obtain the measurement result q - which also plays 
the role of the normalization factor of the post-measurement wave function - reads: 

p{q)= J G,{q-qmq)\'dq. (2) 

As a matter of fact, sharp (direct) von Neumann position measurements are the idealized 
special case while unsharp measurements - though not necessarily with the Gaussian profile 
- are the ones which we encounter in practice and which suite a tractable theory of real-time 
monitoring of the position of a single quantum particle. 

In our discretized model of monitoring a single particle, we assume an unknown initial 
wave function tpolq) and we are performing consecutive unsharp position measurements 
of resolution a at times t = r, 2r, . . . , resp., yielding the corresponding sequence qt of 



measurement outcomes. Between two consecutive unsharp measurements the wave function 
evolves according to its Schrodinger equation (self-dynamics). 

The resolution a and the frequency 1/r of unsharp measurements should be chosen in 
such a way as not to heavily distort the self-dynamics of the particle. It turns out that the 
relevant parameter is a^r, we call 



the strength of position monitoring. If stands for the characteristic extension of the 
wave function, then 70"^ is the average decoherence rate at which our monitoring distorts 
the monitored particle's self-dynamics. We should keep this rate modest compared to the 
rate of the Schrodinger evolution due to the Hamiltonian H of the monitored particle. 
Low values of the strength 7 may, however, result in low efficiency of position monitoring 
and slow convergence of our method of wave function estimation, cp. Sec IIIII The above 
constraints on a^r can in general be matched with further ones - see Sec. IV - that assure 
the applicability of the continuous limit and its analytic equations. 



II. MONITORING THE WAVE FUNCTION 



While it seems plausible that after a sufficiently long time t the sequence of unsharp po- 
sition measurements provides enough data to estimate it may come as surprise that 
position measurements enable a faithful monitoring of the full wave function ipti^) as well. 
Let's just outline the reason. Measuring the position q at times t = r, 2r, ... on an system 
with evolving Schrodinger wave function ipt is equivalent to consecutive measurements of the 
Heisenberg observables qt = ex'p{itH)q exp{—itH) on a system with static wave function tpQ. 
The set of Heisenberg coordinates {qt} will exhaust a sufficiently large space of incompatible 
observables so that their measurements will lead to a faithful determination of tpo and - this 
way - to our faithful determination of tpt for long enough times t. In the degenerate case 
H = 0, monitoring turns out to be trivial: For long enough times, a large number t/r of un- 
sharp position measurements of resolution a is equivalent with a single sharp measurement 
of resolution oj^/Tfr^ position monitoring thus yields just preparation of a static sharply 
localized wave function - an approximate 'eigenstate' of q. 

Our monitoring of the wave function means a real-time estimation of it, where the quality 
of monitoring depends on the fidelity of the estimation. We start from a certain initial 
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estimate ipQ and simulate its evolution according to the self-dynamics of the particle, which 
is assumed to be known, until time t = t. Immediately after we have learned the first 
position Qt- from the first measurement on the particle, we update the estimate according to 
the same rule ([T]) as the actual wave function of the particle and renormalize it: 

ip^{q) — ^ normalization x \/Ga{q — Qt)'^t{'1) ■ (4) 

This update resembles the Bayes principle of non-parametric statistical estimation. We 
repeat this procedure for t = 2t,3t . . . to expect that the estimated V't and the observed 
wave function ipt will converge! A rigorous proof of convergence is missing. In the continuous 
limit, nonetheless, it has been proved for the general case [2^ - not excluding the lack of 
convergence in specific degenerate cases when the set of Heisenberg observables {qt} remains 
too narrow to determine ipo- This is the case for example for a two-dimensional separable 
dynamics in two coordinates x, y, where only the coordinate x is monitored. Rather than 



pursuing the rigorous theoretical conditions of convergence cf. 2]J], we turned to numerical 



tests of continuous measurements that have definitely confirmed our method. 

III. NUMERICAL SIMULATIONS 

We simulated the evolution of a single hydrogen atom subject to continuous measure- 
ments in several potentials. However, the conclusions of our discussion are not restricted to 
hydrogen atoms; similar results can be expected for atoms with higher masses in appropri- 
ately scaled potentials. 

The coupled evolutions of wave function, measurement readout and the estimated wave 
function were simulated numerically by discretizing the corresponding stochastic differential 
equations (cp. Sec JIVp . For this purpose we employed a corresponding scheme of Kloeden 



and Platen which is accurate up to second order in the time step of the discretization [l8l. |22|. 

In order to study the relation between the evolution of the wave function of the particle 
on one hand and the evolution of its estimate on the other hand we first restrict to a one- 
dimensional spatial motion. In this case the graphical representation is the simplest and 
thus gives a clear picture of the convergence between real and estimated wave functions. 
As example we consider a hydrogen atom situated in a quartic double well potential with a 
shape as depicted in Fig{T] 
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FIG. 1: Time sequence of real probability density of atomic position Itpl"^ (blue solid line) and the 
estimated probability density of position (read dashed line). The solid black line represents 

the double well potential as a function of position. Its minima are 189 ^m apart and the height of 
the central maximum is given by 1 x 10~^^ eV. 



We assumed a continuous measurement of the position of the hydrogen atom with strength 
7 = 9.9856/(/im)^s. In order to get an impression of the dimensions of the measurement, 
let us invoke Eq. ([3]) to note that this value of 7 may correspond, e.g., to single Gaussian 
measurements with spatial resolution of a = 1.4mm repeated at time periods r = 50ns. The 
spatial resolution of the single weak measurements 23| is thus 140 times poorer than the 
width cr = 10/im of the initial Gaussian wave function of the atom. In Fig{T]we have depicted 



snapshots recorded at different times of the spatial probability density 



of the H-atom 



(blue solid line) and the squared modulus |'?/^'^(x)p of the estimated wave function (red dashed 
line). Initially both probability densities assume the form of Gaussians which differ in width 
and location within the double-well potential. The sequence of pictures demonstrates the 
convergence of the densities in the course of a continuous measurement. 



The real probability density 



of the H-Atom, which possesses initially a slightly 
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higher mean energy than the middle peak of the potential, oscillates back and forth between 
the sides of the potential. The oscillatory motion of the centre of would also be 

expected qualitatively without measurements -as well as from a classical particle of the same 
mass moving in the double well. However, the real probability density does not spread as it 
would do without measurements. This localisation effect caused by the continuous unsharp 
position measurement adapts the motion of the H-Atom to that of a classical particle which 
is perfectly localised at each instance. This illustrates the influence of the measurements 
and points to a particular kind of control of the wave function that can be exercised by 
means of unsharp position measurements. A smaller measurement strength would lead to 
less disturbance of the unitary motion in the potential but also to less speed of convergence 
between real and estimated probability density. 

The estimated probability density which is centered initially on the right-hand 

side of the potential, follows the real wave function until after approximately one oscillation 
period the corresponding probability densities coincide and evolve identically thereafter. 
But not only the probabilities to find the atom at a certain position converge, in fact the 
complete wave function ip^x) and its estimate iIj^{x) coincide after a sufficiently long period! 
It can be proved analytically, that the estimation fidelity F = \ which measures the 

overlap between the wave functions ip and ip'^, when averaged over many realisations of the 



continuous measurement increases with time and converges to 1 [20|] . Numerical simulations 
for the double- well potential show that in a typical realisation with initial values as described 
above and measurement strength 7 = 9.9856/ (/im)^s the fidelity amounts to more than 95% 
after 1.5 oscillation periods. 

Figj2] shows the evolution of estimation fidelity for of a continuously measured H-atom 
moving in a plane under the infiuence of a mexican-hat potential. The latter is a rotationally 
symmetric version of the double well in two spatial dimensions. For the sake of simplicity 
we assumed that both position coordinates are simultaneously and independently measured 
with the same strength 7. Such a continuous measurement of both coordinates typically 
yields evolutions of the estimation fidelities which are shown in FigE] for several values 
of the measurement strength 7. In all depicted cases the fidelity comes very close to one 
within a period of 5ms, i.e., estimated and real wave function then coincide. Thereafter 
the dynamics of the wave function including the infiuence of the measurement can thus be 
monitored with perfect fidelity. 
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FIG. 2: The fidelity of estimation F is plotted as a function of the duration of the continuous 
measurement of a H-Atom moving in a mexican hat potential for different values of measurement 
strength 7 in units of 10/(/im)^s. The height of the mexican hat's central peak situated at the origin 
of the reference frame is given by 1.07 x lO^^^eV, its minima lie on a concentric circle with radius 
40/xm around the origin. The wave function -0(3;) and its estimate ip^{x) are initially Gaussians 
centered at (— 55^m,— 14.8/xm) and (— 103.6/im,— 103.6/im) with widths of 10/um as well as 5/im, 
respectively. The plots demonstrate that perfect fidelity is reached eventually for all considered 
measurements strengths. However, the convergence time decreases with increasing 7. 

One might doubt that our monitoring remains efficient for heavily complex wave functions 
like those developing in classically chaotic systems. Instead of the integrable Mexican hat 
potential, this time we study the chaotic Henon-Heiles potential which depends on the radius 
r as well as on the azimuth 0: 

y (x, y) = A [r^ + ar^ + br^ cos(30)] . (5) 

For this potential we simulated continuous position measurement and monitoring with the 
following results. The saturation of fidelity is reassuring: the estimate converges to the 
real wave function (Fig|3]), which is found at a just slightly longer time scale than in the 
integrable Mexican hat potential, whereas the wave functions show an apparently irregular 
complex structure. FigJH shows the estimate and the real wave functions with an average 
overlap (i.e. a fidelity) of 91.58%. Indicating already an accurate overall estimation of the 
real wave function, this value still includes in small areas of the potential differences of the 
corresponding probability densities up to 35% of their highest peak. In particular, FigjH 
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FIG. 3: The fidelity of estimation F is plotted as a function of the duration of the continuous 
measurement at strength 7 = 12.351/(/xm)^s of a H-Atom in the non-integrable Henon-Heiles po- 
tential ([5]) with parameters A = 5.44 x 10~^'^eV/ (/im)^, a = 13.09;um^ and b = 36.18/xm. The wave 
function tp{x) and its estimate ^^(x) are initially Gaussians centered at (— 14.8^m,— 29.6/um) and 
(— 29.6/im,— 29.6;Lim) both with widths of 10/im, respectively. We find that the fidelity converges 
to 1 and therefore our estimate becomes a good approximation of the real wave function. 



indicates that faithful monitoring is not only possible when the shape of the wave function 
of the particle is close to a Gaussian but also for rather complex shapes. 

Monitoring, i.e. continuous unsharp observation, has a specific capacity. It is its ro- 
bustness against external unexpected perturbations. To demonstrate such a robustness, we 
assumed that close to saturation of the estimation fidelity, like in FigJH our atom in the 
Henon-Heiles potential is suddenly perturbed, e.g., by a collision with an environmental par- 
ticle (here another hydrogen atom). For simplicity, we assume a momentum kick along 
the x-direction which implies a multiplication of the real wave function by the complex func- 
tion exp{ipxX / h) hence the estimated wave function has to start a new cycle of convergence. 
We map the momentum kick to a temperature by ksT = p^/m as if it had a thermal 
origin, just to give a hint of its strength. Numeric results (Fig|5l) show that the estimation 
fidelity recovers against these momentum perturbations (cp. movie 2J]). In realty, repeated 
random perturbations might prevent perfect monitoring and fidelity will saturate at less 
than 1. This case is beyond the scope of our present work, its study will be of immediate 
interest since real systems are subject to various noises that are not measured at all. The 
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FIG. 4: The real (blue) and estimated (red) spatial probability densities IV'P and IV'^P) whichever 
is the bigger one, are depicted in the Henon-Heiles potential after time 3.15ms, at fidelity 0.9158, 
for the same initial states as in FiglSj 

monitoring theory at non-optimum efficiency has been outlined earlier [2^. 




FIG. 5: At t = 3.15ms and fidelity 0.9158, exactly when the snapshot of FigSlwas taken, the H- 
atom is hit by another thermal H-atom which causes an immediate drop of fidelity. We supposed 
a single momentum transfer px = \/mkBT in the rr— direction only, at different temperatures. The 
fidelity recovers and converges to 1 after some time which depends on the temperature, i.e., on the 
strength of the momentum kick. 
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IV. THE ITO-METHOD 



The discrete sequence of unsharp measurements (fT][2!) and wave function updates 



possess their continuous hmit 25|] if we take r ^ and cr^ ^ cxd at I/7 = rcr^ = const. In 
this 'continuous hmit' both the true wave function ipt{(l) and the estimated wave function 
ipti(l) become continuous stochastic processes such that they are tractable by two stochastic 
differential equations respectively. The position measurement outcomes qt do not yield a 
continuous stochastic process themselves. It is their time-integral Qt, specified below, that 
becomes a continuous stochastic process. 

Let us consider the discrete increment of the true wave function during the period r, cf. 
([T]). In Dirac formalism we get: 

1 



AlV') = exp(-M)-^Va(g-g)|^) - |^) . (6) 

For simplicity, we omit notations of time dependence t. The symbol (q) stands for {ipt\q\'ipt) ■ 
In the continuous limit, the Eq. ([6]) transforms into the following Ito-stochastic differential 
equation 



d\ij) = (^-zH - l{q - {q)yyt\ij) 



+ Y^q-{q)){dQ-{q)dtm (7) 

The equation of the discrete increment A|?/;^) of the estimate (slight change) assumes the 
same form as Eq. ([6]) of but the normalization factor differs from l/p{q), cf. (jlj). Yet, 
it yields the same Ito-stochastic differential equation as the equation above. The estimated 
state l?/;^) must be evolved according to the same non-linear differential equation ([7]) that 
describes the evolution of the monitored particle's state {ipt)- These two equations are 
coupled via the stochastic process Q whose discrete increment is defined by AQ = qr, in 
the continuous limit this means formally Qt = /q qsds where is the measured position 
at time s. In realty, the random process Qt is obtained from the measured data {qt}- 
If the measurement is just simulated, like in our work, then in the continuous limit AQ 
transforms into the Ito-differential dQ whose random evolution can be generated by the 
standard Wiener process W via dQ = {q)dt + 'j~^^'^dW. Of course, Qt breaks the symmetry 
between the stochastic processes ipt and ipt because dQ/dt fluctuates around {'ipt\(l\''Pt) and 
not around 

11 



The stochastic differential equation ([7]) - combined with the same one for is a 

suitable approximation of our discrete model (Sees. II-III) under two conditions: (i) a single 
measurement does not resolve any particular structure of the wave function, i.e., a ^ cr^ 
where is the width of the spatial area on which ip is not negligibly small. Thus a^, can, 
e.g., be of the order of magnitude of the available width of the confining potential, (ii) The 
length of the time period r between two consecutive measurements is small compared to 
the timescale of self-dynamics generated by the Hamiltonian H . Then the discrete model of 
position monitoring and wave function estimation becomes tractable by the time-continuous 
equation ([7]) depending on the single parameter 1/7 = rcr^, cf. also Eq. ([3]). 

V. SUMMARY 

We simulated numerically continuous position measurements carried out on a single quan- 
tum particle in one- and two-dimensional potentials. In order to monitor the evolution of an 
initially unknown state of the particle in a known potential, we estimated its wave function 
and updated the estimate continuously employing the measurement results. 

Our simulations show that for all considered potentials the overlap between estimated 
and real wave function comes close to 1 after a finite period of measurement -guaranteeing 
thereafter precise knowledge of the particle's state and a real-time monitoring of its further 
evolution with high fidelity. The power of our method is indicated by the ability to monitor 
even the motion of a particle in a classically chaotic potential subject to continuous position 
measurement. 

We thus demonstrated, that monitoring the complete state of a quantum system with 
infinite dimensional state space is feasible by continuously measuring a single observable on 
a single copy of the system. Moreover, the simulations indicate that our monitoring method 
is robust against sudden external perturbations such as occasional random momentum kicks. 
How much and what kinds of external noise this monitoring scheme tolerates is important 
for its applicability in control and error correction tasks, and might be object of future 
research. 
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